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INPUT  IMPEDANCE  OF  AN  INFINITE  LINEAR  SCANNED  ARRAY  OF  FLAT  RECTANGULAR 
STRIPS  INSIDE  A  DIELECTRIC  LOADED  PARALLEL  PLATE  WAVEGUIDE 

I.  INTRODUCTION 

In  most  phased  array  studies  the  main  cbjeotive  has  been  to  achieve 
a  wide  angle  impedance  match,  preferably  with  a  wide  bandwidth.  Different 
configurations  and  techniques  have  been  proposed  for  this  purpose.  For 
instance,  Magill  and  Wheeler  [1],  and  also  Munk  with  his  co-workers  [2] 
proposed  different  dielectric  loading  for  achieving  wide-angle  impedance 
match.  The  geometry  of  the  problem  can  play  an  important  role  in  satisfy¬ 
ing  some  requirements.  For  example,  a  parallel  plate  waveguide  can  be 
used  advantageously  to  study  the  H-plane  and  E-plane  scan3  independ¬ 
ently.  With  these  general  objectives  of  wide-angle  and  wideband  H-plane 
performance,  this  paper  is  aimed  at  deriving  an  analytical  expression  for 
the  input  impedance  of  an  infinite  linear  array  of  flat  rectangular  probes 
(monopoles)  inside  a  dielectric  loaded  parallel  plate  waveguide  using  some 
simplifying  assumptions,  so  that  theoretical  results  can  predict  the  cor¬ 
responding  experimental  data  reasonably.  The  parallel  plate  waveguide  is 
short  circuited  at  a  distance  z«  -L  from  the  array  of  probes  [Fig.  1]. 
The  raonopoles  are  flat  rectangular  strips  fed  from  coaxial  lines.  The 
most  difficult  point  of  this  kind  of  analysis  is  in  the  determination  of 
the  current  distribution  along  the  probes.  This  can  be  done,  in  prin¬ 
ciple,  by  solving  an  integral  equation,  which  involves  a  tedious  task  of 
numerical  computations.  In  order  to  alleviate  this  difficulty,  it  has 
become  common  to  assume  a  sinusoidal  current  distribution  along  a  current 

element  [2j-[73.  This,  one  of  the  simplifying  assumptions  used  here,  is 
Manuscript  approved  May  27,  1986. 


an  approximation  and,  therefore,  any  theoretical  analysis  based  on  this 
can  only  approximately  '..redict  the  corresponding  experimental  results. 
Another  simplifying  assumption  consists  of  neglecting  the  effect  of  the 
gap  at  the  junction  between  the  probe  and  its  coaxial  feed  line  [4,5]. 
One  of  the  results  of  this  study  is,  therefore,  to  find  how  much  agreement 
between  these  approximate  theoretical  results  ana  the  corresponding 
experimental  data  is  achieved.  This  will  be  shown  mi  inly  tnrough 
graphical  display. 

With  the  above  objectives  and  assumptions  in  mind,  a  mathematical 
formulation  of  the  problem  which  includes  several  dielectric  sheets,  is 
presented  in  section  II.  In  section  III,  a  sinusoidal  current  distribu¬ 
tion  it  introduced  and  then  the  unknown  complex  amplitudes  of  the  field 
components  are  determined.  A  method  to  improve  the  accuracy  of  thi3  cur¬ 
rent  distribution  by  using  the  stationary  property  of  tne  input  impedance 
is  also  discussed.  A  theoretical  expression  for  the  input  impedance  of 
the  infinite  linear  array  of  the  probes  is  derived  in  section  IV.  In 
section  V,  various  special  cases,  including  the  vor? ^3pondence  between  the 
infinite  linear  array  and  its  equivalent  simulators,  are  discussed. 
Various  theoretical  and  experimental  results  obtained  from  different 
sources  are  compared  in  section  VI.  This  section  also  includes  the  com¬ 
parison  between  theoretical  and  experimental  results  of  the  input  impe¬ 
dance  of  an  infinite  linear  array  of  flat  rectangular  probes  inside  a 
parallel  plate  waveguide  loaded  with  and  without  dielectric  sheets. 
Experiments  were  performed  using  several  waveguide  simulators.  However, 
results  involving  only  one  simulator  are  presented  for  economy  of  space. 
It  was  found  that  the  theory  based  on  the  simplifying  assumptions 
mentioned  earlier  can  predict  the  general  or  qualitative  behavior  of  input 
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impedance  or  VSWR.  It  ha a  also  been  shown  that  a  dielectric  sheet  with  a 
relatively  low  dielectric  constant  can  improve  the  bandwidth  and  match  of 
the  array.  Munk  and  his  coworkers  [2]  had  also  shown  that  a  sheet  of  low 
dielectric  constant  can  improve  scanning  characteristics.  The  relative 
position  between  the  dielectric  slab  and  the  array  chosen  in  ref.  [2]  is, 
however,  different  from  that  considered  here.  The  frequency  range  used 
for  this  investigation  was  about  3*0  GHz  to  M.O  GHz,  over  which  the 
experimental  data  show  a  good  impedance  match.  This  frequency  range 
inc  cates  a  bandwidth  of  28$  with  respect  to  the  center  frequency  3.5  GHz 
having  VSWR  <2.  The  corresponding  theoretical  result  shows  somewhat  lower 
bandwidth. 

II.  MATHEMATICAL  FORMULATION  OF  THE  PROBLEM 

The  geometry  of  the  problem  is  3hown  in  Fig.  1 .  It  consists  of  a 
parallel  plate  waveguide  of  width  b  and  short  circuited  at  a  distance  z= 
-L  from  an  infinite  linear  array  of  y-directed  flat  rectangular  probes  (or 
monopole)  located  in  the  plane  z-0.  In  the  positive  z  direction,  the 
waveguide  extends  to  infinity.  The  probes  are  separated  a  distance  a  in 
the  x-direction.  All  elements  Ci.e.,  probes)  of  the  array  are  excited 
with  equal  amplitude.  The  phase  shift  per  unit  length  is  a  constant  h, 
which  is  a  function  of  frequency  and  scan  angle.  In  this  formulation,  it 
is  assumed  that  the  array  is  inside  a  dielectric  slab  with  relative 
dielectric  constant  e2.  which  again  is  sandwiched  between  two  other  slabs 
having  relative  dielectric  constants  ei  and  eg,  respectively.  These 
dielectric  slabs  are  iueluued  for  better  matching.  The  primary  reason  for 
choosing  flat  rectangular  probe,  instead  of  circular,  is  for  achieving 
better  bandwidth.  Since  the  width  of  the  probe  is  very  small  compared  to 
the  wavelength  Jn  the  medium  in  which  the  probe  is  embedded,  the  current 
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in  the  probe  can  be  assumed  to  he  confined  to  the  y-direction  (along  the 
length  of  the  probe)  only.  The  electromagnetic  field  satisfies  the 
following  time  harmonic  Maxwell’s  equations. 


7  x  -jay  H 
0 


(1) 


7  x  K  -  juc^E  +  yQJy 


(2) 


where  yQ  is  the  unit  vector  in  the  y-direction.  The  relative  dielectric 
constant  is  a  piecewise  constant  having  different  values,  in  general, 
in  different  regions  of  the  z-dimension.  The  assumed  harmonic  time 
dependence  expCjat)  is  suppressed  for  convenience.  The  effect  of  the  gap 
at  the  junction  between  a  probe  and  its  coaxial  feed  is  neglected. 

The  x-dependent  part  of  the  current  on  the  probe  array  can  be 
represented  by  an  infinite  number  of  delta  functions  in  the  following  way. 


-jhna 


5(x-x'-na) 


(2a) 


where  the  factor  exp (-jhna;  represents  the  uniform  phase  shift  in  the  x- 
direction.  Since  the  expression  (2a)  is  a  periodic  function,  it  can  also 
be  represented  by  a  Fourier-series  of  the  following  form. 


I  e  ^hna5(x-x'-na) 


2im, 


I  l  e-J(h^)(x-x') 


n»"» 


(2b) 


The  right-hand  side  of  (2a)  implies  that  the  x-dependence  of  the  fields 
will  also  have  a  similar  infinite  series  representation  containing  a 
factor  exp[-j(h+2nir/a}xj.  It  may  be  noted  that  this  expression  is  also 
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equivalent  to  a  Floquet  mode  [7]  or  an  unnormalized  mode  function 
associated  with  a  one  dimensional  grating  [8]. 


Since  the  electric  current  is  in  the  y-direction,  in  absence  of  any 

dielectric  slab  the  electromagnetic  field  could  have  been  expressed  in 

terms  of  a  single  component  Ay  of  a  vector  potential  A  [4],  However, 

this  approach  is  not  valid  in  the  present  situation  where  dielectric  slabs 

are  present.  Therefore,  a  more  general  method,  namely-  modal  analysis  [8] 

appropriate  for  a  linear  infinite  array  will  be  adopted  here.  In  this 

method,  let  us  introduce  two  transverse  scalar  mode  functions,  d>  (x,y) 

nm 

and  i}>n[n(x,y).  The  word  transverse  means  x.y-dimensions  (which  is  perpen¬ 
dicular  to  the  propagation  direction  z) .  These  functions  are  defined  by 


♦nm'1*’ 


(3a) 


(3b) 


where 


g  (  1  when  m*0  i 
oo  0  when  m^O 


(30 


*  h  +  ?mr/a  (3d) 

Ktnm  ’  tK  J*  <3e) 

Note  that  the  functions  given  by  (3a)  and  (3b)  have  the  same  x-dependence 

shown  in  (2b).  Furthermore,  <j>nm(x,y)  and  ^  (x.y;  vanish  at  y  -  o  and 

oji  nm 

v  »  b.  Define  two  vector  functions  by 


b 


(4a) 


* 

nm 


(x,y) 


-  7 


t  nm 


(x,y) 


*nm(x'y)  “  ~Vnm(x’y)  (4b) 

A  ft  A  ^  A  A 

where  V.  *  y.  r-  +  y^  r-  ,  x^  and  y^  are  unit  vectors  in  the  respective 

t-  0  dX  O  dy  O  O 

coordinate  directions.  Note  that  I  (x,y)  and  $  (x,y)  have  the  follow- 

nm  nm 

ing  properties. 


7t  \  ■  ■  -KtnB2  KmU’y) 

7t  '  t*m(y’x>  x  2oJ  ’  0 

Vt  •  ?n«(x-7)  ' 

't  ’  t?raU’y)  x  ^o1  ‘  0 


(5a) 

(5b) 

(6a) 

(6b) 


where  z  is  the  unit  vector  in  the  z-direction. 
o 


II 


♦nm(x,y)  •  $* ._'(x,y)  dxdy  -  &  ,6  , 

nm  n  m  nn  mm 


;; 


7nm(x,y)  •  Vm'(x,y)  dxdy  -  6  ,6  - 

nm  n  m  nn  mm 


// 


Vx*y)  *  zo x  *5'm'(x»y)<!xdy  - 0 


(7a) 

(7b) 

(7c) 


The  asterisk  signs  in  (7a)  to  (7c)  indicate  respective  complex  con¬ 
jugates.  The  symbols  6  .  and  6  ,  represent  Kronecker  delta  functions 

nn  mm 

with  the  usual  meaning.  The  region  of  integration  in  (7a)  to  (7c)  is  the 
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3  3 

unit  cell,  i.e.,  -  and  o  S  y  S  b.  With  these  orthonormal  vector 
mode  functions  $nm(x,y,)  and  ^nm^x,yi^’  the  comPiete  electromagnetic 
field  can  be  represented  in  the  following  manner  [8]: 


Et(r) 


•  l  l 

n— •  o-o 


V™(z)5  (x,y)  - 

nm  nm 


l 

n— • 


m-o 


V*E(z)$  (x,y)xz  (8a) 
nm  nm  o 


l  I™(z)z  X*  (x,y)  +  l  l  I^(z)l  (X,y)  (8b) 

nm  o  nm  _  nm  nm 


Jueo£iEz(r) 


l 

n-  -• 


m-o 


iS(2)kLw*-i,) 


(8c) 


ja1u0Hz(r) 


l 

n-  -• 


l 

m-o 


vJE(z ) K?  ♦  ( x ,y ) 

nm  tnra  run 


(3d) 


The  subscripts  t  and  z  denote  transverse  (i.e.,  x,y)  and  z  components  of 

TM  tm 

the  respective  field  vectors.  The  functions  V  (z)  and  I  (z)  are  the 

nm  nm 

TE 

modal  voltage  and  current  amplitudes  for  TM  waves.  Similarly,  V  (z) 

nm 

TE 

and  I  (z)  belong  to  the  TE-type  waves.  Adopting  the  procedures  similar 
nm 

to  that  in  ref.  [8],  it  can  be  shown  that  these  modal  voltage  and  current 
amplitudes  obey  the  following  transmission  line  differential  equations. 


-  ~  v™(z)  -  r.  z™ 
dz  nm  inm  inm 


tTM,  . 
I„ m(z) 
no 


(9a) 


.TM, 


~  l“‘(z)  -  r  Y™  V™(z)  +  ff  y  •  #*  (x,y)J  (r)dxdy 
dz  nm  inm  inm  nm  ''  o  nm  y 
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(9b) 
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(10a) 


where 


VTE(z)  -  jwu  ITE(z)  -  r  ZTZ  ITE(z) 
dz  nm  J  o  nm  inm  inra  nm 


I™(2)  -  y  (10b) 


,TE  ,.TE 


dz  nm  inm  inm  inm  '  J  J"o 


^  ■  t‘L  <^'n 


(11a) 


K0  -  »  ✓  Uoe0 


(11b) 


i  -  1,2,3  and  H 


(11c) 


~TE  .  /viE  . 

z.  *  i/Y.  *  Juju  /r. 

inm  inm  J  o  inm 


modal  impedance  for  TE-wavea 


(lid) 


-TM  , /VTM  „  ,, 

Z.  *■  1/Y.  -  T.  /jut  e, 

inra  inm  inm  J  o  i 


-  modai  impedance  for  TM-wave3 


(lie) 


Jy(?)  -  -~fl)  6(z)  -  Jy(y)6(z) 


(12) 


for  |x|  i  and  o  iy  id 


W  -  width  of  the  flat  rectangular  probe  in  the  x  direction 
d  -  effective  length  of  the  probe. 
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Note  that  the  probe  is  infinitely  thin  in  the  z  dimension.  The  integra¬ 
tions  in  (9b)  and  10b)  are  over  the  surface  of  the  probe  as  defined  in 
(12).  The  propagation  constant  Tinm  given  by  (11a)  is  different  in 
different  regions  and  when  it  is  real,  the  corresponding  waves  are 

evanescent.  For  propagating  modes  (such  as  main  beam  and  grating  lobes) 

TE  TM 

T.  is  purely  imaginary.  For  convenience,  let  us  introduce  U  and  U 
inm  mn  nm 

defined  respectively  by 


jTE 

nm 


>u0JJ> 


»*  (x.y) 
nm 


Z(y) 

W 


dxdy 


(13a) 


and 


J™ 

nm 


2nm 

jVoe2 


J/i 


** m<x»y> 

nm 


¥l± 

w 


dxdy 


(13b) 


In  order  to  solve  for  V™(z),  I™(z),  vJE(z),  and  lIE(z), 

nm  nra  nm  nm 


it  i3  con- 
,TM, 


venient  to  obtain  a  second  order  differential  equation  for  each  V  (z)  and 

nm 

TE  TM  TP 

V*  (z).  This  can  be  d  me  by  eliminating  I  (z)  and  I  (z)  from  (9a), 
nm  nm  nm 

(9b)  and  (10a),  (10b).  The  results  are  (using  (12),  (13a)  and  (13b) 


and 


-  -r  c)  V™(z)  -  -Ulfl6(z) 
2  inm  nm  nm 


,TM. 


dz 


(14a) 


dz 


-r  *)  VTE ( z )  -  -  U1!i6(z) 
2  inm  nm  nm 


.TE. 


(14b) 
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These  two  equations  must  be  solved  subject  to  the  following  boundary 


I 


condi v Ions. 

At  the  source  (i.e.~  at  z-o) 

V™(z*o*)  «  V™(z-o") 


(15a) 


(z-c-0 


(z-o  ) 


nm 


(15b) 


(z«o+) 


(z-0  ) 


(16a) 


(z»o+) 


»2(*-°’) 

nm 


(16b) 


TM  TM  TP  TP 

Since  there  i3  no  otiie-  sources,  V  (z),  I  rz),  V  (z)  and  I*  (z)  are 

nm  nm  nm  nm 

TM  TP 

continuous  at  z  «  6  and  5.  ♦  5,«  I  (z)  and  I  (z)  can  be  deter- 

c.  s  s  3  nm  nm 

TM  xr 

mined  from  V  '  (z)  and  V  ~(z),  respectively,  via  Eqs.  (9a)  and  (10a).  In 
nm  nm 

addition, V*^(z)  and  V*^(z)  must  vanish  at  z  -  -L.  For  z>6  +6,,  VTS(z) 
nm  nm  s  3  nm 

TM 

and  V  ‘  (z)  behave  like  outgoing  waves. 
nm 

The  solutions  of  (l^a)  and  (l^b)  subject  to  the  boundary  conditions 
discussed  above  can  now  be  expressed  in  the  following  manner. 
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A™  3inh{r,  (s+L) }  -UzS-i0  (17a) 

nas  ■  i  nm  2 

3™m[e-p!r2,,m2!~  4K(™,‘e*p,‘r2M(2V*,,l'"VzSo  (17b) 

V™(2)  ’  \  3™m!exp!‘r2nn2i*"ifr(r’>n!'expl'r2„B<263-z)i'  <17o) 

C™!exp('r3ns2>*  S3!«(n'“’'e,!p,'r3nm(23s*2{3‘2,|!’  SsSzSVS3 

(17d) 

TM 

D‘”exp(-rj._mz)  z£5  (1?e) 

na  ^nm  3  3 


A™3inh  |ri™tz*ut 


-LSzS'£2  (18a) 


8ILtel!p(r2n»l)'4frCn'B,'exp|"r2iim(2V2,l)’  ~i  SzSo 

(18b) 


^nm(z)  “  S  B2Ltexp(-W)+4ff(n'ffl),exp{'*r2nm(2V2)'1-1' 


Cl  8c ) 


4mteXp(’r3nmz)+R3,4(n'ra),exp{-r3nm(26S+253-z)^'  VzSV*3 


TE 

Dnmexp(-r«nm2) 


3ntn  s  3  '  * "  '3 - 3 

(l8d) 


(i8e) 
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Rgff(n,m)-[R]1>2(n,3!)+exp{-2rlRm(L-l2)f]/[t+Rj1j2(n,ia).exp{-2rinm(L-l2}] 


H9a) 


R;'i2(n.m).(Elr2nB-e2r,nm)/(Elr2m*c2r,nn)  09b) 

RgJf(n^m)“[R21j3(n.ra)+R^l4(t.,m)‘exp(-2r2n{n53)]/[1+R213(n,m)»R”j,{n?m)‘ 

•  exp  (“2f3nm53)j  (19c) 


^3(n,a)  -  (V3n.-»3r2»)/(V30.  *  £3r2nJ  "9<1) 


93.4(n’")  •  (c3r«r»  -r3nm!/(13r9nm  *  r3n»! 


(I9e) 


B™m  -  C‘+Rllf(n.®)-exp(-2r2nm63)]/[l+R”f(n,ia).R^f(n,m)- 

run 


!X»  ^WW*! 


(19f) 


3Sm  3  BInm *  ^ 1  f { n » m ^  * exP ( -^r.,  £_ )  3 / [  1  +r!.  ( n , m ) • exp ( -?r9  5  ) ]  (19g) 


2nm  Inm  eff' 


“2nm  2 


2nm  a J 


#™  ■  BL-^!OT-exl>t-r2r,m12)t’*neff<n'")!]/teir2ra,eoshIrinm!l'-i‘2i1 


if,  .:n 


(19h) 


Cn”  “  3™m‘exp{'65{r2nm“r3nm)}*[l+Reff(n’m)i/t1"R3!^(n»E)<exp{'2r3nIn65^ 


(19i) 
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°Z  ■  ('♦»3;il(n.=.)]-«pi(Js.S3!(r!lra-r3ran)ll-CnB 


TM 


(19J) 


-  [ftj-  2(n,m)  +  exp  »-2r-^(L-£2J}]/[l+a|  2(n,m)* 


exp  {-2r1nffl(L-£2)}]  (20a) 


aj  ^(n.a)  -  (r1nffi  ~r2ninMr1n2,  *  r2Rm) 


1 


(20b) 


B^ff (nrm)  ’•  [ai  3(n»m)+aJ'  ^nia^exp^ar^o^)  ^‘r|-  ^(n.m)* 


♦exp(~2F_  «„)]  (20c) 

anm  3 


a4  ,(n,m)  •  (r-  -r-  }/(r„  +  r_  ) 

2,3  2nm  3nnr  v  c tun  3nnr 


(20d) 


RM'n,ffl)  *  ^r3nffi“I’J{nm^/^r3na?  +  r%m; 


(20e) 


ts  yiC* 

*  — —  [14f-exp(^r_6„}]/[l^L(n.m)-R^<n,a). 


inra 


m 


2mn  s'J  L  eff' 


eff 


2nra 


(20f) 


B2ta'  B^a-[!'^rf(n>m)'el",('2r’™£»)!/iURi<'<-fn.")-«pC-2r,™«JJ 


2nai  2' 1  1  '  eff 


'2nm  sJ 


(20g) 
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kll  ■  C2n“  •  exP t -r2nml2 )  E '  Cn  ,m )  ] / [ r , nm003h { r ,  nmc L-ic2 ]  n 


(20h) 


(20i ) 


Dm-  [l*4,,(n,»)]-exp|{S3*S3)(r3l)B-r%m)|.C^ 


(20j) 


III.  PROSE  CURRENT  DISTRIBUTION 

The  amplitude  factors  UTE  and  U™  given  by  (13a)  and  ( 1 3b )  which 

nm  nm 

depend  on  the  current  distribution  on  the  probe  still  remain  unknown.  In 

reality,  it  is  difficult  to  calculate  the  exact  current  distribution  in  an 

environment  shown  in  Fig.  1.  Therefore,  various  approximate  methods  are 

sought.  For  example,  one  may  assume  a  known  current  distribution  [3,**. 

5,6]  which  has  been  found  to  be  valid  in  some  situations.  For  a  probe 

this  current  distribution  is  taken  to  be  sin[K(d-y)],  where  d  is  the 

length  of  the  probe,  y  is  the  coordinate  along  which  the  current  varies 

and  K  is  the  free  space  wave  number  of  the  surrounding  medium.  If  b  is 

the  height  of  the  waveguide  in  the  y-direction,  it  has  also  been  shown 

experimentally  [5,6]  that  theoretical  results  based  on  the  assumption  of 

this  sinusoidal  current  are  in  good  agreement  with  those  of  experimental 

for  dS*6b  and  dSl/4,  where  1=2tt/K.  On  the  other  hand,  if  the  probe  is  at 

the  interface  between  two  media  with  propagation  constants  K-j  and  Kg,  then 

the  effective  ourrent  distribution  may  be  taken  [9]  as  sin[Ke(d-y ) ] ,  where 

K  -  /„2  7  Jl  ,  assuming  that  the  permeability  is  the  same  for  both  media, 
e  +  ^2 

All  these  results  are,  of  course,  derived  for  an  antenna  in  an  otherwise 


free  apace. 


Another  approach  to  derive  an  approximate  form  of  current  distribu¬ 
tion  follows  from  the  stationary  property  of  the  input  inpedance  of  the 
probe  with  respect  to  the  current  distribution.  It  is  shown  [4,103  that 
the  input  impedance  of  a  probe  can  be  cast  into  a  variational  form  with 
respect  to  the  variation  of  current  distribution  in  the  probe.  Then  the 
current  distribution  may  be  expressed  in  terms  of  two  or  more  arbitrary 
parameters  which  are  computed  by  requiring  that  the  partial  derivative  of 
the  input  impedance  with  respect  to  each  of  these  unknown  parameters 
vanishes.  When  the  expression  for  the  input  impedance  is  too  complex, 
this  procedure  becomes  tedious. 

Most  of  the  theoretical  models  deal  with  idealized  situations.  For 
example,  for  a  probe,  one  may  assume  a  current  distribution  either  as  a 
known  function  or  a  trial  function  having  a  few  arbitrary  parameters. 
Even  an  integral  equation  for  the  current  may  be  set  up  using  an  idealized 
probe  geometry.  However,  in  actual  practice,  the  probe  geometry  compris¬ 
ing  the  coaxial  feed  system  differs  from  the  idealized  theoretical  config¬ 
uration.  As  a  result,  the  actual  current  may  depart  substantially  from 
the  theoretically  computed  one.  Therefore,  in  order  to  minimize  the  com¬ 
plexity  in  theoretical  computation,  the  following  approach  will  be  chosen 
here. 

The  current  distribution  on  the  probe  will  be  assumed  to  be  known 
and  has  the  following  form  (see  Eq.  12). 

I(y)  “  Iosin[Ke(d-y)]  (21) 
oSySd,  |x|$W/2 
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where  Kg  is  the  effective  propagation  constant.  If  the  medium  in  which 

the  probe  is  embedded  has  a  relative  dielectric  constant  e2,  then  Kg  = 

io  -  K0/~e^  .  On  the  other  hand,  if  the  probe  lies  at  the 

interface  between  two  media  with  relative  dielectric  constants  e.  and  e., 

i  J 

then  K  =  K  /  (e, +e.)/2  . 
e  o  i  j 

TE  TM 

Using  now  (21),  the  amplitude  factors  U  "  and  U  can  be  calculated 

nm  nm 

from  Eqs.  (13a)  and  (13b).  The  results  are: 


„  K  .  2-6  „  K  S  F 

„TE  „  .  2(  ed-u  omi1/2f  xn  >  n  era 

U  «*  -21  up  K  sin  l-s—Jl — r—J  It: - j  “ - 

nm  o  o  e  v  2  ab  '  VK.  ' 

tnm  K2 
em 


( 22a) 


,TM  2IoKe  „  2  ,  2fKedu2  >1/2  fmir/b>  SnFem 

'»  " '  3=vJ  •3in  ‘T%) 


(22b) 


F  =  1  -sin2(nnrd/2b)/sin2(K  d/2) 
em  e 


(23a) 


Sn  ‘  3in  'Kxn"/2>/K„“/2> 


(23b) 


(mTT/b)2  -Ke2 


(23c) 


IV.  INPUT  IMPEDANCE  CF  THE  INFINITE  LINEAR  ARRAY  OF  PROBES 

The  input  impedance  Zin  is  defined  [3,^3  by  the  ratio  of  the  complex 
power  radiated  by  the  probe  and  the  mean  square  input  current,  i.e., 
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mum  mmm 


-  j  jj  St(x,y,7.=o)«y  ^  ,)v(y)ds/[^|  J(o)  j2]  .  (24) 


The  integration  is  over  the  probe  surface  inside  the  unit  cell 

a 

containing  the  origin.  Jy(y)  is  given  by  02)  and  (?1).  The  electric 
field  Ej/x.y.z-o)  can  be  calculated  from  (8s)  using  (17c-).  0 8c) ,  (19g), 
(20g),  (22a),  and  (22b).  Carrying  cut  the  indicated  operations  in  (24), 
the  input  impedance  Z^n  can  be  expressed  as 


OT  00 


Zin  *  TST5"  tan2CKed/2)  l  l 

n»- «  m**o 


(2-6  )/k  \2/f  S  \2 

v  om;  I  xn  1 1  em  n  1 

r2nm  '  tnm*  \  K2  / 
em 


A 

nm 


JZ  K 


ab  tar^(V/2)  I  I  Cf^)2  r2ml 


|OsJ-Vt" 


(24a) 


where 


T~  - 
nm 


[  1  +R£ff  (n,m)  -exp(-2r5mn6 J ][ 1  -P^^n.rn)  -expC-Zr^iJ  ] 


2nm  s 


eff' 


[ 1  +4ff<n,m)  R^rf.(n,m) ,exp{-2r9rim(i^*6a) } ] 


2nm  2; 


eff 


2nm  2  sJ 


(24b) 


r11- 

nm 


[UR^r(n,m)-exp(-2r2mi53)][l-R^f(n,m)-exp(-2r2nmil2)] 


:n 


rRef f (n ’ m 5  * Ref f ( n ’ ra ; ‘ exp ^“2r2nm U2+6s ^ ^ 


‘11 


(24c) 


Z 

o 


/  uQ/e 


o 


(24d) 


17 


WSWWMPjSIW? 


Phased  arrays  are  generally  designed  in  such  a  way  that  only  the 
lowest  order  mode  or  the  main  beam  (for  which  n*-o=m  in  (24a))  pro¬ 
pagates.  However,  due  to  dielectric  loading  the  first  few  higher  order 

modes  (i.e.,  grating  lobes,  such  as  when  n-±1 ,  m=o;  n**o,  m»1 ,  etc.)  may 
become  propagating  waves  inside  a  dielectric  layer.  Even  in  this  case, 

Finm  is  real  for  large  value3  of  n  and  m  (which  represent  evanescent 

I  1  * 

modes).'  Therefore,  for  large  values  of  n  and  m,  both  T-1-  and  T  tend  to 

nm  nm 

approach  unity. 

From  Eq3.  (3d),  '3e),  (Ha).-  and  23a)  to  (23c),  one  fines  that  Kxn, 

Ktnm  and  r.  behave  like  n  as  n  tends  to  infinity  for  a  fixed  m.  On  the 

other  hand,  for  a  finite  value  of  W,  Sn  Is  on  the  order  of  1/n  as 

When  n  remains  fixed  and  m  becomes  very  large  K\  and  r,  tend  to  m. 

tnm  inm 

A 

Fem  and  K*n  approach  unity  and  m2 ,  respectively,  as  m  assumes  a  large 
value.  From  this  observation,  it  follows  that  the  first  tern  (double 

3  *7 

series)  of  (24a)  converges  like  l/n-’and  1/m'as  n  and  m  tena  to  infinity, 
respectively.  The  second  term  of  (24a),  however,  converges  like  \/xr>  and 
1/m^  as  n  and  m  approach  infinity,  respectively.  On  the  otherhand,  if 
(KxnW)  is  very  small,  the  impedance  series  converges  slowly  with  respect 
to  n. 

V.  SPECIAL  CASES 

A.  IN  ABSENCE  OF  ANY  DIELECTRIC  SLAB 
In  this  case  »  1 


r.  -  r  «  /  „  2  v  2 

inm  nm  K.  -  K 
tnm  c 

Reff(n,m)  *  *>.«)  *  0 


1 ,2,3,4 
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~  I 

R~ff (n,ra) 


ll  -  x11 

nm  nra 


Rglf(n,m)  -  exp{-2rnm(L-i2)} 

1  -exp  (-21*  L) 
nra 


in 


jZ  K  K 
J  o  o  e 

2ab 


-  (2-<$  )  K 

tan^Uyi/2)  l  l  ■  — 

n=*-°  m*o  1  nm  tnm 


)2(4ELH)a[1-exp(-2rnmL)] 


“5f 

K2 


em 


-jZ  K 
_  o  o 

ab 


tan“(K  d/2) 
e 


n*-a 


I 

m-o 


rrair/b 

LK„ 

tnm 


)2r 


nm 


(4£i)2[i-exP(-2r  l)] 


K2 

em 


nm 


(25) 


In  this  case,  K0  may  be  taken  as  K  -  u>  /  yem  .  In  the  absence  of  any 

— 2r  L 

short  circuit  at  z  *  ~L,  one  may  take  e  nm  «  o.  This  is  Justified  by 
taking  the  limit  L  with  the  assumption  that  Rernm>c  for  all  n  and 
m.  For  propagating  modes,  this  assumption  implies  Inclusion  of  slight 
loss  in  the  medium. 

B.  PROBE  AT  THE  CENTER  OF  THE  BROAD  FACE  OF  A  RECTANGULAR  WAVEGUIDE 

In  Fig.  1 ,  a  is  taken  to  be  the  distance  between  the  consecutive 
probes  of  a  linear  array  in  the  x-dimension.  This  is  also  shown  in  Fig. 
2a  for  convenience.  The  width  of  the  parallel  plate  waveguide  is  b  in  the 
y-dimension.  In  this  special  case,  the  rectangular  waveguide  has  dimen¬ 
sions  a  and  b  in  the  x  and  y  directions,  respectively. 

The  center  of  the  probe  is  at  x-a/2.  This  probe  at  the  center  of  a 
rectangular  waveguide  is  equivalent  to  a  linear  array  of  probes  whose 
adjacent  elements  are  180°  out  of  phase  in  a  parallel  plate  waveguide  as 
shown  in  Fig.  2b.  This  observation  follows  from  image  theory.  Therefore, 
when  the  uniform  phase  shift  h  associated  with  the  probes  in  Fig.  1  or 
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Fig.  2a,  is  given  the  value  ir/a,  the  problem  of  Fig.  1  becomes 
equivalent  to  that  in  a  rectangular  waveguide  with  a  probe  at  the  center, 
x  *=  |  .  As  a  result,  the  following  transformations  take  place. 

f?xn  »  h  +  2mr/a  -  (2n+1  )ir/a  (26a) 


K.  »  [k2  *  (mir/b)2]1/*-[{2n+1  )ir/a}2+(mir/b)2]1/2  (26b) 

unm  Aii 

rinm  *  [  l(2n+1  )ir/a}2+(mir/a)2-K2ei]  /  .  (26c) 

It  can  readily  be  seen  from  (24a)  to  (24c)  and  (26a)  to  (26c)  that 
the  argument  of  the  summation  with  respect  to  n  associated  with  the  input 

impedance  Zin  becomes  then  a  function  of  (2n+1)2.  If  it  is  assumed  that  G 

2 

is  a  function  of  (2n+1)c,  then  it  can  be  shown  easily  that 

00  00 

l  G[(2n+1)2]  -  2  l  G[ (2n+1 )2]  .  (27) 

n— "  n«o 

Using  this  information,  the  desired  input  impedance  of  a  centrally 
located  rectangular  probe  in  a  rectangular  waveguide  is  given  by 
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'in 


1Z  K  K  d  -  «  •  (2-6  )  /  K  F  S  \ 

0  °  e  tan2(K  d/2)  l  l  V-22  I  J2L2y!_  \ 

dW  6  n*°  m"°  r2nm  V  Klm  / 

'  tnm  sin  / 


1*1 

nm 


J2Z  K  2  tan2 (K  d/2)  »  »  /  (mir/b)F  S 

oe  e  rfl  eni  r 

"X“ 


ab  K  e0 

O  2 


®  »  /  (mir/bjF  S  \ 

I  I  ( - ^  ) 

n=0  \  Ktnm  Kem  > 


r  t11 

2nm  nm 


(28) 


In  (28)  the  relations  (26a)  to  (26c)  are  to  be  used. 

The  relation  (28)  can  further  be  specialized  for  the  situation  when 
there  is  no  dielectric  3labs  inside  the  rectangular  waveguide,  i.e.,  when 
c^me2me^m  e^-1 .  Following  the  procedure  described  in  Section  V-A,  we 
have  the  desired  expression  obtained  from  (28). 


JZ  K 


•i.  ^  tan2(K0d/2)  j_  J  V[l-exp(-2rnBL] 

\  tnm  em  ' 


n-o  m«o  nm 


J2ZCK0  . _ Zr 


ab 


.  •  •  /  (mir/b)F  S  \ 

tan2(K  d/2)  l  I  (  ;  emn  )  rnm[l-exp(-2rn|]1L)] 

n-o  m-1  \  K.  K*  / 

'  tnm  em  '  ^9) 


The  Eq.  (29)  can  be  simplified  to  the  following  form. 


-1Z  K  _  -  •  (2-«  )  /F  S 

Ztn  •  ^  «o4/2)  1  1 

n»o  m-o  ....  . 

em 


-«  )  /F  S  \ 

[i-xp(-2r„mi.)] 

™  \  K».  / 


(30) 


It  may  be  noted  that  situations  represented  by  Eqs.  (28)  and  (29)  or 
(30)  correspond  to  a  particular  scan  angle  for  a  given  frequency. 
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C.  SOME  SIMULATOR  MODELS 


Since  It  is  very  expensive  and  time  consuming  to  build  a  periodic 

array  with  a  large  number  of  elements  for  experimental  testing  and 

measurements,  use  of  simulator  models,  which  are  simple  in  structure  and 
contain  only  a  few  elements,  has  been  proposed  and  proved  to  be  successful 
[11,12].  The  principle  involved  in  this  simulation  is  also  discussed  ir 
these  references  [11,12]  and,  therefore,  will  not  be  repeated  here.  One 
of  the  restrictions  in  simulator  application  is  that  each  simulator 
represents  only  one  sample  of  scan  conditions  in  the  actual  array. 
Consequently,  several  simulators  are  generally  built  for  a  study  of  a 
given  array.  For  example,  the  special  case  discussed  in  section  V-5 

consisting  of  a  rectangular  waveguide  with  a  probe  at  the  center, 

represents  a  particular  scan  condition  defined  by  h**Kosin8*ir/a  for  the 
array  problem  depicted  in  Fig.  1.  The  H-plane  scan  condition  for  this 
infinite  linear  array  can  be  simulated  approximately  by  the  simple  TE10 
mode  in  a  rectangular  waveguide.  The  number  of  elements  N,  to  be  placed 
in  the  rectangular  waveguide  and  the  broad  dimension  a  of  this  guj.de 
depend  again  on  the  scan  angle  9,  satisfying  the  condition 


h  -  K  sin9  *  Tf/a 
o 

or  sin8  =*  A /2a 


(31a) 


If  dx  is  the  interelement  spacing,  then  the  number  of  elements  required 
for  a  given  scan  angle  is  approximately  given  by 


N  -  a/d 

x 


A 

2dx  sin8 


(31b) 
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For  a  given  scan  angle,  say  0-35°,  the  required  broad  dimension  of  the 
rectangular  waveguide  is  a-0.8717A.  If  dv  is  chosen  as  d  -0.575A,  then 
N=ta/dx'*1 .516.  Thus,  the  approximate  number  of  elements  is  1.5,  i.e.,  one 
element  is  bisected  by  one  wall  of  the  guide  (in  the  x-direction)  and  the 
other  element  is  placed  at  a  distance  dx*2a/3  from  that  bisecting  wall. 
Since  the  bisected  element  is  short  circuited  by  the  wall,  it  is  observed 
that  its  existence  (or  non-existence)  has  no  effect.  This  observation 
then  implies  that  a  rectangular  waveguide  with  a  probe  at  a  distance  2a/3 
from  one  wall,  a*0.877A  being  the  broad  dimension  of  the  guide,  when 
excited  by  the  TE10  mode,  simulates  the  scan  condition  (0-35°)  of  an 
infinite  array  of  probes  in  a  parallel  plate  waveguide.  For  making  this 
observation  more  instructive,  the  input  impedance  of  an  off-centered  probe 
in  a  rectangular  waveguide  is  computed  (see  Appendix)  separately  and  the 
corresponding  numerical  results  are  compared  with  the  appropriate  situa¬ 
tions  for  an  infinite  array  probe,  for  which  Eq.  (25)  has  been  used.  The 
equivalence  between  any  other  simulator  and  the  infinite  linear  array  in  a 
parallel  plate  waveguide  for  a  given  scan  condition  can  be  explained 
similarly  [9,10]. 

VI.  DISCUSSIONS  OF  NUMERICAL  AND  EXPERIMENTAL  RESULTS 

Over  the  past  several  decades  to  the  present,  numerous  theoretical 
and  experimental  investigations  of  phased  arrays  composed  of  electric 
dipoles  or  monopoles  have  been  made.  In  theoretical  analysis  a  knowledge 
of  the  current  distribution  along  these  dipoles  or  monopoles  is 
essential.  However,  it  is  not  known  a  priori.  It  can  be  determined  in 
principle  by  solving  an  integral  equation  for  the  current.  Once  the 
current  distribution  is  determined,  other  parameters  of  interest  (such  as 
input  impedance,  VSWR,  radiation  pattern,  etc.)  can  be  computed  from  the 
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knowledge  of  the  current.  The  solution  of  the  integral  equation  involving 
current  is  generally  accomplished  numerically.  Since  this  procedure  is 
tedious  and  time  consuming,  one  generally  assumes  a  sinusoidal  current 
distribution  in  theoretical  analyses  for  simplicity.  Although  such  an 
assumption  simplifies  the  analysis  considerably,  differences  are  expected 
to  exist  between  the  theoretical  and  the  corresponding  experimental 
results.  Nevertheless,  this  simplified  assumption  may  be  considered 
satisfactory,  provided  the  theoretical  results  can  predict  at  least  the 
qualitative  behavior  of  the  input  impedance  and  other  observables  of 
practical  interest,  so  that  the  experimental  results  can  be  better 
understood.  This  i3  one  of  the  objectives  of  thi3  theoretical  study.  It 
may  be  noted  that  for  some  problems  [5,6]  assumption  of  known  sinusoidal 
current  distribution  produced  theoretical  results  which  were  in  good 
agreement  with  those  of  the  corresponding  experiments. 

Before  we  proceed  to  make  a  comparison  between  various  theoretical 
and  experimental  results,  let  us  digress  for  a  moment  in  order  to  present 
the  interesting  result  to  be  discussed  in  the  next  section. 

VI.  A.  CONNECTION  BETWEEN  CIRCULAR  AND  FLAT  RECTANGULAR  PROBE  IMPEDANCE 

OR  VSWR 

One  may  choose  either  a  circular  probe  or  a  flat  rectangular 
probe  for  a  phased  array  composed  of  monopoles.  A  question  then  arises 
whether  there  is  a  simple  connection  between  two  corresponding  observables 
of  the  same  kind  (such  as  input  impedance  of  the  probe,  VSWR,  etc.).  This 
information  is  particularly  important  when  some  results  (theoretical  or 
experimental)  using  one  kind  of  probe  are  known,  then  how  the  theoretical 
or  experimental,  as  the  case  may  be,  parameters  for  the  other  kind  of  the 
probe  may  be  chosen  so  that  the  two  corresponding  results  should  be  the 
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3ame.  It  appears  that  h.A.  Wheeler*  was  also  interested  in  this  matter 
and  derived  a  relation  using  conformal  mapping  of  fields  in  an  unpublished 
report.  Wheeler's  relation  is  W«4p  where  W  is  the  width  of  an  infinitely 
thin  flat  rectangular  probe  and  p  is  the  radius  of  the  circular  probe, 
other  parameters  of  the  two  configurations  (one  with  the  flat  rectangular 
probe  and  the  other  with  the  circular  probe)  remaining  the  same.  The 
implication  of  this  relation  may  be  interpreted  by  using  the  following 
example.  Suppose  there  are  two  identical  rectangular  waveguides,  one  with 
a  circular  probe  of  radius  p  and  the  other  with  a  flat  rectangular  probe 
of  width  W-4P.  These  probes  are  of  the  same  length,  excited  in  the  same 
manner  and  situated  in  an  identical  manner  in  their  respective  wave¬ 
guides.  Then  the  input  impedance,  VSWR,  etc.,  of  the  two  situations  are 
expected  to  be  the  same.  Since  Wheeler's  relation  is  based  on  conformal 
mapping,  valid  for  electrostatic  fields,  it  may  not  be  satisfactory  In  the 
microwave  region.  This  observation  turns  out  to  be  true.  Although 
Wheeler's  relation  is  very  simple,  in  a  frequency  dependent  situation  no 
such  closed  form  relation  can  be  established.  Therefore,  one  shall  have 
to  resort  to  numerical  computations.  In  order  to  obtain  such  a  relation 
let  us  consider  a  known  result  [4]  of  a  circular  probe,  as  shown  in  Fig. 
3,  at  the  center  of  the  broadface  of  a  rectangular  waveguide.  The  wave¬ 
guide  is  short  circuited  at  a  distance  z«  -L  from  the  axis  of  the 
probe.  When  only  the  fundamental  TE1C)  mode  is  propagating,  the  input 
impedance  of  the  probe  is  given  by  [4]  the  following  expression.  It  may 
be  recalled  that  in  this  case  also  the  current  distribution  on  the  probe 
is  assumed  to  be  known  and  is  the  same  as  in  eq.  (21). 


*  This  information  was  communicated  to  the  author  by  Dr.  W.  K.  Kahn. 
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of  the  second  kind,  of  order  zero.  In  the  absence  of  the  short  oircuit 
(i.e.,  L+«)  the  corresponding  expression  for  the  input  impedance  of  the 
probe  car.  be  shown  to  have  the  following  fo^rc. 
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The  result  corresponding  to  eq.(32)  for  a  flat  rectangular  probe  is  given 

by  (30).  In  the  absence  of  any  short  circuit  the  input  impedance  can  be 

obtained  from  (30)  by  letting  exp(-2r  L)  -  o.  Numerical  values  of  these 

nm 

impedance  expressions  are  obtained  using  the  parameters,  a  *  6.0  em,  b  « 
3.4036  cm,  L  *  2.21234  cm  and  ®,  d  -  1.651  cm,  1.75  cm  and  1.90  cm.  The 
radius  p  of  the  circular  probe  is  kept  fixed  at  0.127  cm.  However,  the 
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width  W  of  che  flat  rectangular  probe  is  varied  from  Up  to  5p.  The  results 
are  shown  in  Fig.  3a  thru  4c  in  the  form  of  VSWR  as  a  function  of  frequen¬ 
cy.  The  voltage  standing  wave  ratio  (VSWR)  is  computed  from  the  relation 


■ 

VSWR  =  1 


V50  -1 1 


Z.  /50+1  i 
in 


(34) 


It  is  found  from  these  numerical  results  that  the  input  impedance  and  VSWR 
of  a  circular  probe  remain  practically  unchanged  when  it  is  replaced  by  a 
flat  rectangular  one,  provided  that  W  -  4,882p,  ether  parameters  being 
the  3ame.  This  relation  should  not  be  considered  exact.  It  may  change 
with  frequency,  length  of  the  probe,  distance  L  of  the  short  circuit, 
etc.  These  observations  are  based  on  the  numerical  results  shown  In  Fig. 
3a  thru  4c.  However,  the  relation  W  -  4.882p  appears  to  be  adequate  for 
the  frequency  range  and  other  parameters  considered  here.  Fcr  the  purpose 
of  comparison,  the  VSWR  corresponding  to  Wheels.'  's  relation  W  -  4p  is 
also  presented.  Although  Wheeler’s  relation  is  obtained  by  using  con¬ 
formal  mapping,  valid  for  electrostatic  fields,  results  (such  as  impe¬ 
dance,  VSWR,  etc.)  obtained  by  using  his  relation  approach  closer  to  those 
obtained  by  using  the  relation  W  *  4.882p,  when  the  frequency  is  higher 
and  the  probe  length  is  longer. 

The  relationship  W  =  4.882p  implies  that  the  expressions  for  the 
Input  impedances  of  flat  rectangular  probes  developed  in  this  paper  for 
various  situations  can  also  be  used  lor  the  corresponding  cases  where  the 
former  probe  is  replaced  by  a  circular  one  of  radius  W/4,882.  It  is,  of 
course,  assumed  that  other  parameters  remain  unchanged. 
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VI.  8.  A  CIRCULAR  PROBE  IN  AN  INFINITELY  LONG  RECTANGULAR  WAVEGUIDE 
WITHOUT  SHORT  CIRCUIT 

Before  presenting  our  experimental  and  theoretical  results  for  the 
linear  array  of  flat-  rectangular  monopolea  in  a  parallel  plate  waveguide 
or  its  equivalent  simulator,  let  us  first  make  a  comparison  between  some 
experimental  and  theoretical  results  related  to  Collin’s  [4]  and  Al- 
Hakkak's  [5]  model,  since  their  various  assumptions  are  similar  to  ours. 
This  model  consists  of  a  circular  probe,  fed  from  a  coaxial  waveguide  in 
an  infinitely  long  rectangular  waveguide  without  short  circuit.  The  main 
reasons  for  choosing  this  model  here  are  the  following: 

1.  The  theoretical  results  can  easily  be  derived  from  Collin’s  work 
[4]  with  a  miner  modification.  Thi3  was  done  by  Al-Hakkak  [5], 
who  also  performed  experiments  and  compared  with  the  theory. 

2.  The  theoretical  results  for  this  model  can  also  be  obtained 
approximately  from  our  theory  for  a  simulator  using  the  relation 
W  =  4.882p  in  eq.  A-10  in  view  of  the  finding  in  the  previous 
section. 

Although  Chang  and  Kahn  [6]  studied  a  similar  problem  using  a  flat 
rectangular  stub  inside  an  infinitely  long  rectangular  waveguide  without 
using  any  coaxial  feed  line,  their  experimental  and  theoretical  results 
cannot  be  used  here  directly  for  comparison.  This  model  does  not  have  any 
problem  of  gap  at  the  coax  junction. 

The  theoretical  expression  for  the  input  impedance  is  given  in  eq. 
(33),  which  is  a  slight  modification  of  the  result  given  by  Collin  [4], 
In  view  of  the  relation  discussed  in  section  VI-A,  the  numerical  values  of 
eq.  (33)  for  a  given  probe  radius  p  can  also  be  obtained  approximately 
from  eq.  ( A— 1 0 )  with  W  *  4.882p,  and  the  omission  of  the  exponential 
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factor  in  (A-10),  other  parameters  remaining  the  same.  For  the  use  of  eq. 
(33),  it  should  be  remembered  that  only  the  fundamental  TE10  mode  is 
propagating.  In  addition,  the  theory  does  not  take  into  account  the  gap 
(or  aperture)  at  the  junction  of  the  probe  and  the  coaxial  connector  and 
assumes  a  known  sinusoidal  current  distribution  along  the  probe.  A 
further  limitation  of  the  theory  is  the  assumption  that  the  diameter  of 
the  probe  is  very  small  compared  to  the  wavelength.  In  ref.  [5],  Al- 
Hakkak  finds  also  that  the  length  of  the  probe  d  should  be  such  that 
d<0.6b  and  d<A/4,  where  b  is  the  shorter  dimension  of  the  broadface  of  the 
waveguide.  Chang  and  Khan  [6]  also  discussed  these  limitations  of  the 
theory.  It  may  be  noted  again  that  the  sinusoidal  current  distribution 
assumed  here  and  in  refs.  [4,5]  is  also  used  by  Stark  [3]  among  others  for 
similar  theoretical  study  of  phased  array  antennas  having  monopole  or 
dipoles  as  elements. 

Figures  (5a),  (5b)  and  (5c)  show  3ome  theoretical  and  experimental 
values  of  normalized  (with  respect  to  50Q)  input  impedance  of  circular 
probes  having  three  different  diameters,  but  of  the  same  length  in  an 
infinitely  long  rectangular  waveguide  without  short  circuit.  Although  the 
parameter  c,  the  outer  diameter  of  the  coaxial  line,  does  not  enter  into 
the  theory  due  to  the  approximation  made,  it  does  affect  experimental 
results.  The  probes  are  located  at  the  center  of  the  broadface  of  the 
waveguide.  The  results  show  a  general  agreement  between  the  theory  and 
the  experiment.  The  main  reasons  for  the  existence  of  difference  between 
the  theoretical  and  experimental  results  may  be  attributed  to  the  gap  (or 
aperture  effect)  and  the  assumption  of  a  sinusoidal  current  distribution 
along  the  probe.  In  order  to  observe  the  effect  of  the  gap  alone,  the 
probe  was  removed  from  the  coaxial  connector  in  an  experiment.  It  was 
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then  noticed  that  there  was  hardly  any  transmission  of  energy  through  the 
gap  into  the  waveguide  (showing  extremely  high  VSWR) .  This  phenomenon 
might  have  been  known  also  to  others  and,  perhaps,  motivated  them  to 
neglect  the  effect  of  the  gap  in  similar  theoretical  works  as  a  good 
approximation.  This  may  not  be  justified  since,  when  the  gap  and  the 
probe  are  present  simultaneously,  the  field  distribution  at  the  junction 
is  altered  due  to  the  coupling  between  the  gap  and  the  probe.  In 
addition,  the  combined  effect  of  the  gap  and  the  probe  diameters  appears 
to  play  a  role.  For  example,  an  increase  in  probe  diameter  causes  an 
increase  in  the  input  impedance  and  thereby  improves  VSWR.  An  increase  in 
the  probe  diameter  is  also  associated  with  a  decrease  in  the  resonance 
frequency  (at  which  the  input  impedance  is  real)  of  the  system.  This  was 
also  observed  by  Al-Hakkak  [5].  It  is  also  observed  (not  presented  here) 
both  theoretically  and  experimentally  that  an  increase  in  probe  length  d 
improves  the  VSWR  within  a  limited  frequency  band.  It  may  be  noted  also 
that  at  frequencies  above  4.22  GHz,  the  chosen  probe  length  d  -  1.778  cm 
does  not  satisfy  the  condition  d<A/4. 

In  spite  of  all  these  restrictions  made  in  the  theory,  Al-Hakkak's 
experimental  results  are  surprisingly  in  good  agreement  with  the  theory. 
Figure  6d  shows  theoretical  results  corresponding  to  the  parameters  used 
by  Al-Hakkak  [5J,  These  are  essentially  Al-Hakkak’s  results  (recomputed). 
The  agreement  is  so  close  that  the  experimental  and  theoretical  results 
are  not  shown  separately  in  Fig.  5d.  It  may  be  mentioned  here  that  Chang 
and  Khan  [6]  also  obtained  good  agreement  between  their  theory  and 
experiments  in  which  there  is  no  gap  effect.  Their  theory  is  also  similar 
to  that  of  Collin’s  where  use  of  known  sinusoidal  current  distribution  is 
made.  Note  that  Chang  and  Khan  [6]  apparently  determined  the  actual 
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current  distribution  by  solving  an  integral  equation  numerically «  They 
found  that  the  actual  current  distribution,  so  determined,  does  not  differ 
much  from  a  sinusoidal  current  distribution  provided  d<A/4. 

The  results  of  this  section  may  be  considered  as  precursors  of  what 

can  be  expected  of  the  results  for  the  linear  array  of  monopoles  in  a 

parallel  plate  waveguide  and  hence  may  be  used  for  guidance.  The 

experiment'll  results,  Figs.  5a  to  5c,  were  generated  in  our  laboratory. 

VI.  C.  INPUT  IMPEDANCE  OF  AN  INFINITE  LINEAR  ARRAY  OF  FLAT  RECTANGULAR 
MONOPOLES  INSIDE  A  PARALLEL  PLATE  WAVEGUIDE 

In  this  case,  eqs.  (24a)  and  A-8)  are  used  for  theoretical 
computations.  Equation  (24a)  is  valid  for  both  the  presence  and  absence 
of  a  dielectric  slab  and  for  all  scan  angles.  On  the  other  hand,  eq.  (A- 
8)  is  used  for  the  case  without  any  dielectric  and  for  those  scan  angles 
and  frequencies  3hown  in  Fig.  6a.  In  the  presence  of  a  dielectric  slab 
only,  the  eq.  (24a)  is  used  with  eim1*E2“e4  £j**1.3»  a-2.84« (2.54)cm-3d1 , 
where  2.84 * (2. 54) cm  is  the  longer  dimension  of  thr  broadface  of  the 
simulator  waveguide,  d1  is  the  offset  distance  of  the  center  of  the  probe 
from  one  end  of  the  simulator  waveguide  [Fig.  A-1],  b-1 .34*(2.54)cm, 
W»0.2’(2.54)cm,  d-0.702’ (2.54)cm,  L-0.863* (2.54)cm,  6  -0.984* (2.54)cm, 
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6^*1 .0*(2.54)cm.  All  the  measurements  were  taken  using  a  simulator  (Fig. 
A-1 )  with  3d.|=*its  longer  dimension  of  the  broadface.  Note  that  a  of  Fig. 
1  is  not  the  same  as  a  of  Fig.  A-1.  In  the  absence  of  any  dielectric, 
eqs.  (24a)  and  (A-8)  provide  identical  results,  subject  to  the  condition 
d1 -2.84* (2.54)  cm/3-distance  between  the  consecutive  probes  of  the  array 
in  the  parallel  plate  waveguide.  This  is,  of  course,  expected. 

Figure  6b  shows  the  theoretical  and  experimental  results  in  the 
absence  of  any  dielectric  slab.  Here  also,  the  theoretical  and  experi¬ 
mental  results  are  in  qualitative  agreement.  Similar  results  with  a 
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dielectric  slab  are  shown  in  Fig.  6c.  The  experimental  results 
show  a  tuning  effect  in  the  presence  of  the  dielectric  slab.  Although 
this  effect  is  not  pronounced  in  the  theoretical  result;  it  indicates  such 
an  effect  qualitatively.  It  may  be  noted  that  the  experimental  results 
show  better  bandwidth  and  match  than  the  corresponding  theoretical 
resuits.  Some  theoretical  computations  (not  presented  here)  also  show 
that  the  bandwidth  and  .r-ntch  cun  be  improved  to  some  extent  by  decreasing 
the  inter-element  s  <r.oing  a.  The  reasons  for  the  differ*'-”®  between  the 
theoretical  and  experimental  results  are  the  same  as  tu  .e  discussed  in 
the  previous  v.  ■  VI. B.  It  shou.’  1  now  be  apparent  that  discussions  and 
presentation  of  theoretical  and  experimental  results  in  section  VI. B 
provide  better  understanding  of  the  results  of  this  section. 

VII.  CONCLUSIONS 

For  theoretical  analysis  of  problems  consisting  of  electric  dipole 
or  raonopole  elements,  it  is  common  practice  to  assume  a  sinusoidal  current 
distribution  along  such  current  elements.  This  is  done  mainly  to  simplify 
the  complexity  of  the  problem.  Although  such  an  assumption  is  inaccurate, 
it  can  predict  qualitatively  expected  behaviors  of  various  experimental 
results  (input  impedance,  VSWR,  etc.).  These  expectations  were  Illus¬ 
trated  by  comparing  theoretical  and  experimental  resO.ts  obtained  here  and 
from  previously  published  literature.  It  was  one  of  the  alms  of  the 
present  study  to  make  such  comparisons.  The  primary  causes  of  the 
difference  between  the  theoretical  and  experimental  results  are  explained. 

In  addition,  certain  experimental  results  show  a  better  bandwidth 
(about  28?)  and  match  than  the  corresponding  simplified  theoretical 
results.  It  is  also  shown  both  theoretically  and  experimentally,  that  a 
slab  of  low  dielectric  constant  (such  as  1.3,  I.1*)  prope'-ly  placed  in 


frcnt  of  the  array,  can  improve  the  bandwidth  and  VSWR  at  wide  scan 
angles.  Munk  and  his  coworkers  [2]  had  also  shown  that  a  sheet  of  low 
dielectric  constant  can  improve  scanning  characteristics.  The  relative 
position  between  the  dielectric  slab  and  the  array  chosen  in  ref.  [2]  is, 
however,  different  from  that  considered  here.  The  frequency  range  used 
for  this  study  is  about  3.0  GHz  to  4.0  GHz. 

Another  interesting  result  is  the  relation  between  the  width  W  of  a 
flat  rectangular  probe  and  the  radius  p  of  a  circular  probe  situated 
similarly  inside  two  identical  rectangular  waveguides.  The  probe 
impedances  and  VSWR's  for  these  two  situations  are  found  to  be  the  same 
for  all  practical  purposes  when  W  is  properly  related  to  p.  The  proper 
relation  found  here  is  W**4.882p  for  a  given  set  of  waveguide  parameters 
and  frequencies.  This  factor  was  obtained  numerically  representing  a 
correction  to  the  electrostatic  result  W»4p  found  by  Wheeler. 

It  may  be  observed  from  this  study  insofar  as  the  simplified  theory 
is  concerned,  that  it  provides  results  relatively  easily,  and  does 
correctly  indicate  the  trend  of  the  experimental  results. 
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comparison  here. 
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APPENDIX 


RADIATION  FROM  AN  OFF-CENTERED  PROBE  INSIDE 
A  RECTANGULAR  WAVEGUIDE  SHORT  CIRCUITED  AT  A  DISTANCE  L  FROM  THE  PROBE 

The  geometry  of  the  situation  is  shown  in  Fig.  A-1.  It  consists  of 

a  rectangular  waveguide  of  dimension  a  and  b  (aSb)  and  short  circuited  at 

a  distance  z«  -L  from  an  off-centered  flat  rectangular  probe  of  length  d 

and  width  w.  The  axis  of  the  y-directed  probe  lies  at  x-d.,  and  z»o.  In 

the  positive  z-direction,  the  waveguide  extends  to  infinity.  For 

simplicity,  it  is  assumed  that  there  are  no  dielectric  slabs  inside  the 

waveguide.  The  current  distribution  is  the  same  as  given  by  (21). 

Although  one  can  use  a  method  similar  to  that  presented  in  section  2,  the 

absence  of  any  dielectric  permits  one  to  use  only  a  single  component  of 

the  vector  potential  A.  Since  the  current  is  directed  along  y,  one  finds 
-►  * 

A**y0Ay.  Then  it  can  be  shown  that 

E  *  -juA  +  W*A/(jujeouo)  (A-la) 


H  =•  VxA/u 

0 


(A-lb) 


Therefore, 


E 

y 


(A-2a) 


(A-2b) 
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Let  us  construct  a  Green's  function  G(r,r')  which  satisfies  the  following 
differential  equation  and  boundary  conditions  [4]. 

(7  +Ko2)G(rtr')  -  5(x-x')5(yv~y')<S(z) 

r  »(x ,y,z)  and  r'  *  (x',y',o)  (A-3) 


G  -  o  at  x«o,a  and  z— L 
3G/3y  *  o  at  y-o,b 


(A-4) 


G  also  satisfies  radiation  condition  at  z»®. 

Then 

A  -  -u  J/G(r,r')j  (y')ds  (A-5) 

y  °  s  3 

The  integration  is  over  the  surface  of  the  probe  at  z»o,  and  Jg(y)  *• 
I(y)/W,  where  I(y)  is  given  by  (21).  The  Green's  function  G(r,r')  can 
be  expressed  C1)]  in  the  fol^uwing  form 


G(r ,r ' )«  - 


1 


00 

ab  ^ . 
n-l  m*o 


«  (2-6  ) 

r  otn '  ,  rnirx->  ,  rmrx'i 

L  — ? -  sin  — j  sin  [— — J- 


nm 


•cos  (2-^jcos  (^-)  [  exp  { -fnm  |  z  | } -exp  { -rnm ( 2L+z ) }  ] 


where 


rnm'  [(nir/a)2+(mir/b)2-K02]  1 


(A-6) 
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In  this  case  (see  Eq.  (2*0)  the  input  impedance  Zin  is  given  by 


Z.  »  -1/2  jj  E  (x,y,o)  J*(y')ds/[l/2|j  (o)|2]  (A-7) 

in  y  s  s 

3 

The  field  Ey  can  be  calculated  using  (A-2a),  (A-5)  and  (A-6).  Then 
carrying  out  the  operation  indicated  in  (A-7),  the  expression  for  the 
input  impedance  of  an  off -centered  flat  rectangular  probe  in  a  rectangular 
waveguide  can  be  written  in  the  following  form. 


Z. 

in 


%  tan 


K  d 

2  (4-) 


l  l  - 

n-1  m»o 


(2-6  ) 
v  onr 

r 

nm 


F,  S,  nird, 

sin  ( - 1) 


.  [  1  -exp  { -2rnn]L }  ]  (A-8) 


where 


-  t  -[sin  (g*)/sln  t^)]2 


(A-9a) 


*  ,  (-mrw>  .fn'srw-\ 

S1n  *  3in  [2T)/[2T] 


(A-9b) 


K2  -  (ai?/b)2  -K  2 
m  o 


(A-9c) 
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When  the  probe  is  at  the  center  of  the  broad  dimension  of  the  waveguide, 
i.e.,  when  d.j*a/2,  the  expression  for  Z^n  becomes 


in 


jZ  K  tan2(x  d/2)  -  •  (2-6  ) 

°°  o -  II  -ySS-  [l-exp{-2Y  L}]. 

n-o  m=o  nm 


ab 


A  A  , 

*[f.  s.  /k  ]‘ 

1  1m  In  mJ 


(A-10) 


where 


Ynm  *  [  f(2n+1  )ir/a }  ^  *  Kn2]1/2  (A-lla) 

a 

S1n  «  sln{(2n-l)*/a}/{(2n+l)ir/af  (A-11b) 

Note  that  the  expression  (A-10)  is  identical  to  the  relation  (30)  which 
was  obtained  as  a  special  case  of  an  infinite  linear  array  of  probes 
inside  a  parallel  plate  waveguide.  Furthermore,  when  d1  »a/3  or  d-j**2  a/3, 
then  (A-3)  becomes  identical  to  (25)  for  some  specific  scan  angle  and 
frequency  as  mentioned  in  Section  Y-C.  This  will  be  demonstrated  via 
numerical  results.  It  may  be  reminded  that  neglect  of  the  gap-effect  and 
the  assumption  of  known  current  distribution  are  che  usual  approximations 
made  here. 
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PARALLEL  PLATE  WAVEGUIDE  [SIDE  VIEW] 

y 


♦ 


Fig.  1  —  Linear  array  of  flat  rectangular  probes  in  a  parallel  plate  waveguide  for  H-plane  scanning. 
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-3a  -2s  -a  0  a  2a  3a  *9 

Fig.  2a  —  Schematic  arrangement  of  probes  as  in  Fig.  1. 
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Fig.  2b  —  Probes  having  adjacent  phases  180°  apart 
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VSWR 


FREQUENCY  IN  GHz 

Fig.  3a  —  Connection  between  circular  and  flat  rectangular  probes  in  a 
rectangular  waveguide. 

a  -  2.362  •  (154)  cm,  b  -  1.34  •  (2.54)  cm, 

L  *»  0.871  •  (2.54)  cm,  p  0.05  •  (2.54)  cm, 
and  d  0.65  •  (2.54)  cm. 
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THIS  SOLID  CURVE  IS  FOR  A  RECTANGULnn  PROBE 
OF  WIDTH  W=4  x  p= 0.508  cm  [ELECTRO-STATIC  APPROX] 


^  CIRCULAR  I 
PROBEAT  | 
THE  CENTER  1 


THIS  DOTTED  CURVE  IS  FOR  A 
/  RECT.  PROBE  WITH 
W**4.882  x  p=0.62  cm 


THIS  CURVE  IS 
FOR  A  CIRCULAR  PROBE 
OF  RADIUS  P» 0.127  cm 
[.  ROM  COLLIN'S  MODEL] 
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FREQUENCY  IN  GHz 

Fig.  3b  —  Connection  between  circular  and  flat  rectangular  probes  in  a 
rectangular  waveguide. 

a  -  2.362  •  (2.54)  cm,  b  -  1.34  •  (2.54)  cm, 

L  -  0.871  •  (2.54)  cm,  p  -  0.05  •  (2.54)  cm, 
and  d  ™  0.689  •  (2.54)  cm. 
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.  ^  RECTANGULAR  PROBE  OF  WIDTH  w=4  x  p».508  cm 
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FREQUENCY  IN  GHz 


Fig.  3c  —  Connection  between  circular  and  flat  rectangular  probes  in 
a  rectangular  waveguide. 

a  —  2.362  •  (2.54)  cm,  b  «■  1.34  •  (2.54)  cm, 

L  -  0.871  •  (2.54)  cm,  p  -  0.05  •  (2.54)  cm, 
and  d  -  0.748  •  (2.54)  cm. 


43 


FOR  A  RECTANGULAR  PROBE 
OF  WIDTH  w=4  x  p=0.508  cm 


J _ I _ I _ i _ L _ I _ L  1.  ...  1 _ 1  L  I _ L 

2.8  2.9  3.0  3.1  3.2  3.3  3.4  3.5  3.6  3.7  3.8  3.9  4.0  4.1 

FREQUENCY  IN  GHz 

Fig.  4a  —  Connection  between  circular  and  flat  rectangular  probes  in  a 
rectangular  waveguide. 


a  -  2.362  •  (2.54)  cm,  b  -  1.34  •  (2.54)  cm, 
L  —  oo,  p  —  0.05  •  (2.54)  cm, 
and  d  *■  0.65  •  (2.54)  cm. 


4.2 


VSWR 


Fig.  4b  —  Connection  between  circular  and  flat  rectangular  probes  in  a 
rectangular  waveguide. 

a  -  2.362  •  (2.54)  cm,  6  -  1.34  •  (2.54)  cm, 

L  —  oo,  p  «  o.05  •  (2.54)  cm, 
and  d  m  0.689  •  (2.54)  cm. 
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FREQUENCY  IN  GHZ 

Fig.  4c  —  Connection  between  circular  and  flat  rectangular  probes  in  a 
rectangular  waveguide. 

c  -  2.362  •  (2.54)  cm,  b  -  1.34  •  (2.54)  cm, 

L  ™  oo,  p  «■  0.05  «  (2.54)  cm, 
and  d  —  0.748  •  (2.54)  cm. 
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Fig.  Sa  —  Impedance  of  a  circular  probe  in  an  infinitely  long  rectangular 

waveguide. 

a  -  2.84  •  (2.54)  cm,  6  —  1.34  *  (2.54)  cm, 
d  -  0.70  •  (2.54)  cm,  2 p  -  0.049  •  (2.54)  cm, 
and  c  *•  0.163  •  (2.54)  cm. 
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Fig.  5b  —  Impedance  of  a  circular  probe  in  an  infinitely  long  rectangular 

waveguide. 

a  -  2.84  •  (2.54)  cm,  b  -  1.34  •  (2.54)  cm, 
d  «■  0.70  •  (2.54)  cm,  2 p  »  0.125  •  (2.54)  cm, 
and  c  *  0.360  *  (2.54)  cm. 
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Fig.  5c  —  Impedance  of  a  circular  probe  in  an  infinitely  long  rectangular 

waveguide. 

a  -  2.84  •  (2.54)  cm,  b  -  1.34  •  (2.54)  cm, 
d  -  0.70  •  (2.54)  cm,  2 p  -  0.20  •  (2.54)  cm, 
and  c  *■  0.456  •  (2.54)  cm. 
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Fig.  5d  —  Impedance  of  a  circular  probe  in  an  infinitely  long  rectangular 

waveguide. 

a  -  4.76  cm,  b  -  2.22  cm, 
d  “  1.15  cm,  2 p  *■  0.23  cm, 
and  c  —  1.10  cm. 
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Fig.  6a  —  Si -an  anglj  vs  Frequency  [for  a  simulator]. 


a  -  2.84  •  (2.54)  cm, 
b  -  1.34  •  (2.54)  cm, 
dx  -  a/3. 
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Fig.  6b  —  Array  impedance  in  S-band  simulator— Offset  flat  probes 
without  dielectric  sheet 

«i  »  «j  -  «j  -  «4  -  1,  a  ”  2.84  •  (2.54)  cm  “  3</t, 
b\  «  1.34  •  (2.54)  cm,  w  —  0.2  •  (2.54)  cm, 
d  -  0.702  •  (2.54)  cm,  and  Z.  -  0.863  •  (2.54)  cm. 
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Fig.  6c  —  Array  impedance  in  S*band  simulator— Offset  fiat  probe  with 

dielectric  sheet 

*i  **  <2  “  *•»  *  U  *3  “  1.3,  a  *■  2.84  •  (2.54)  cm  **  Zd\, 
b\  «*  1.34  •  (2.54)  cm,  w  »  0.2  •  (2,54)  cm,  d  *  0.702  •  (2.54)  cm, 
L  -  0.863  •  (2.54)  cm,  S,  *  0.984  •  (2.54)  cm.  S3  *  1.0  -  (2.54)  cm 


S3 


